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In this work simple and effective quantization procedure of classical dynamical systems is 
proposed and illustrated by a number of examples. The procedure is based entirely on differential 
' equations which describe time evolution of systems. 

»S 1 Introduction 

The aim of this work is to propose and illustrate by a number of examples a simple method 
of quantization of classical finite-dimensional dynamical systems. Recall that a quantization in 
q-i a broad sense is a transition from classical description of a dynamical system to its quantum 

description. In the early days of quantum mechanics the quantization could be resumed as selecting 
some solutions of classical dynamical equations according to Bohr- Sommerf eld rules. Subsequently 
quantization was understood in a new fashion as replacing of classical observables by operators 
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ading in Hilbert space of states. The prescription for such replacement was formulated in the late 

twenties of the last century and was called canonical quantization. 

According to the canonical quantization procedure pQ, [2] classical coordinates and conjugated 

momenta are replaced by Hermitian operators q l ,pi which subject to the following canonical 

ON . commutation relations 

O 

[q\q j ] = 0, \pi, Pj ]=0, [q\p j \=\h5)l i,j = l,...n. (1) 

X 

Here 8j is the Kronecker delta and I is the identity operator which, as a rule, will be omitted in what 

follows. Other classical observables f(q z ,Pj) are replaced by quantum observables / = f(q\pj)- 
Further, in the context of the procedure quantum dynamical equations originate from the classical 
Hamiltonian equations 

^ = { q \H}, ^ = { Pt ,H}, i = l,...n (2) 

which are written in Poisson bracket form. Recall that the canonical brackets of two arbitrary 
functions / and g are defined by the following formula 



df dg__ dg df 

dpi dq l dpi 



(3) 
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The replacement of the brackets ([3]) by the commutator 



{f,g} — " fi[f>si\ 



(4) 



in the classical equations of motion (j2J) gives the quantum Heisenberg equations 



d<f 
~dt 



1 



dpi 
dt 



1 



i = 



1, ... n 



(5) 



where H = H(q l ,pj) is a quantum Hamiltonian. 

The canonical quantization algorithm works perfectly for most Hamiltonian systems, in par- 
ticular for the systems which Hamiltonian function has the form 



However, for dynamical systems with Hamiltonians more complicated than ([6]), for multi-Hamilto- 
nian systems, for systems which not allows global Hamiltonian structure, and in some other cases 
the canonical quantization often gives indeterminate results or it is non-applicable at all. In 
addition there is evident dependence of canonical quantization results on coordinate choice. 

These reasons have stimulated on the one hand, revision of the canonical quantization algo- 
rithm (for example, geometric quantization [3]) and on the other hand search new quantization 
methods (path integral quantization [3], deformation quantization [5], etc.). At present there are 
a broad spectrum of various quantization schemes and a rich variety of them with the related 
references are listed in the recent review j6]. 

Almost all quantization methods imply the existence of globally defined symplectic, Poisson or 
more rich structures on a phase manifold of a classical system to be quantized. However, certain 
of the dynamical systems does not allow such a structures or these structures are complicated 
enough (for example, dissipative systems). In this case a quantization method which is relied 
solely on classical dynamical equations looks very attractive. Remark, that possibility to deduce 
commutation relations among coordinates and momentum from equations of motion was discussed 
long ago in the papers [7], [8]. Since that time a few works were made in this direction. Here we 
mention the paper [5] where the Weil quantization scheme was employed to obtain commutation 
relations and quantum dynamical equations for some dissipative systems. 

In the present work a simple quantization algorithm is proposed in which classical equations 
of motion play crucial role. Namely, for any classical system the algorithm allows to determine 
simultaneously quantum commutators for initially commutative dynamical variables and quantum 
corrections to the original classical dynamical equations. As a result we obtain commutation 
relations for any pair of basis quantum variables x, y, ... z and equations of motion in the form 




(6) 




k=i 




dx 



k=i 



2 



In the above relations h is Plank constant, g k and r are functions of specially ordered quantum 
variables. These functions can be defined from the conditions that commutation relations are time 
independent and quantum variables are Hermitian. The same conditions allow to obtain quantum 
corrections for integrals of motion if they exist. 

Thus, for any classical system the proposed quantization algorithm gives in our disposal an 
algebra of quantum observables with a set of commutation relations among its base elements 
and a system of quantum evolution equations. To investigate the obtained quantum system it is 
necessary to choose appropriate representation of the quantum observables algebra. This choice 
essentially depends on form of commutation relations and is beyond the algorithm. So in the 
present work we concern only formal algebraic aspects of the quantization procedure and almost 
entirely ignore other ones. 

In the next section the quantization algorithm will be presented in some details and then 
illustrated by several well known examples (charged spinless particle in a magnetic field, anisotropic 
Euler top, etc.). Some remarks concerning the algorithm will be given in conclusion. 



2 Quantization algorithm 

Consider a finite - dimensional dynamical system which evolution is defined by a system of ordinary 
differential equations or, in geometric terms, by some vector field X (see, for example, [in])- More 
precisely, let M be a real n - dimensional manifold and T{M) the algebra of infinitely differentiate 
complex - valued functions on M. From the physical point of view M is a phase manifold and 
J-(M) is an algebra of observables of this system. Time evolution of any function / from J-"(M) 
is determined by the equation 

^ = *(/), f = f(x)EF(M), xeM (7) 

where X is the vector field under consideration. Recall that any vector field is a derivation of 
!F(M). This means that the linearity condition 



X{af + bg)=aX{f) + bX(g) 

and Leibnitz rule 



X(fg) = X(f)g + fX(g) 

hold for arbitrary functions /, g G J-(M) and complex numbers a, b G C. In local coordinates 
x k , k = 1, n, the field X has the form 

n 

X = Y J X\x)d k , d k = —, X k (x) = X(x k )eF(M). (8) 

k=l 

Using the equation (j7j) we can obtain evolution equations for basis dynamical variables which we 
identify here with the coordinate functions x k = x k (x) 



3 



Equations of the form (Q usually derived from experimental data. 

Solution of equations ([HD, at least when the functions X k (x) are analytical, can be represented 
in the form 

x*(t) = x%xl...,x n )= Yl 4 1 ,...,^(4) fcl ---K) fc " (10) 

fcl, ... ,k n >0 

where Xq, i = 1, ...,n are the values of the coordinate functions at t — and s l kl ^(t) are real 
functions of t. This solution defines a one-parameter group of diffeomorphisms of M 

g l x : M -> M, ^x(O) = x(t), ac G M. 

Observables which are invariant under diffeomorphisms from this group are called integrals of 
motion of a dynamical system and can be obtained from the equation 

^M = X(/) = 0, I{x)eT{M). (11) 

To realize the quantization procedure we firstly perform more or less formal transition from 
the algebra of classical observables J-(M) to a commutative algebra J-(M). The transition is 
defined by the map k : x l \— > x l which associates basis elements x l of the algebra J-{M) with 
the coordinate functions x % = x l (x) G JF(M). This map can be extended naturally to the whole 
algebra T{M). Namely, for any function f(x l , ...,x n ) G T{M) we obtain 

K(f(x\ x n )) = /Wx 1 ), K(x n )) = f(x\ x n ) G T{M). (12) 

The above formula will be quite defined if we add to the set of basis elements M = (x 1 , ...,x n ) a 
multiplicative unit I. So let us put 

k(o) = ale T{M), Va G C. (13) 

Now introduce an involution in JF(M) taking into account reality of coordinate functions x % = x l 
where bar denote a complex conjugation in JF(M). Set for any / G T{M) 

<f) = («(/))*• (14) 

Here star denote involution in .F(M). From this formula and the reality x l it follows that the 
basis elements of T{M) are Hermitian 

( x k y = x k . (15) 

Thus J-{M) is a commutative unital C-algebra with involution and its elements are formal series 
in basis elements M = ( 

To represent dynamic of our system in terms of x l define the vector field X counterpart X as 
a derivation in the algebra J-{M) which acts on the basis elements x k according to the formulas 
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X{£ k ) = X(x k ), k = l,...,n, X{I) = 0. (16) 
This relation and the fact that the series 

x\t)=x t (t,x 1 ,...,x n ) = s k..., k Mxl) kl ---(xo) kn i i = h-,n, (17) 

k\ , . . . ,kn 

is a formal solution of the equations 

dx k 



X(x k ) 



dt 

if the series (fTUI) is a solution of the equations ([9]) suggest that the equation (fT8l) can be accepted 
as a dynamical equations for basis elements x k . 

Thus, as a result of our preliminary formal step we obtain (associative) commutative involutive 
C - algebra JF(M) with Hermitian basis elements x 1 and fixed derivation X of this algebra which 
defines the evolution equations (fl8l) . In addition if the original classical system allows integrals of 
motion Ij, j = 1, m < n, we obtain a set of such integrals defined by the equation 

Xfc) = 0, j = 1, m < n, Ij E F{M). (19) 

Now we are ready for the final step of the quantization algorithm, namely for a transition from 
the commutative algebra J~(M) to a noncommutative algebra of quantum observables T{M). Let 
under the transition each basis element x % G T{M) converts in a basis element x % 6 T{M) and the 
unit I converts in its counterpart I. Further, define an order of multipliers in any product of the 
elements x l . Fix, for example, for any monomial in basis elements M = (x 1 , ... ,x n ) the following 
ordered form 

:M Nl ... Nn : = (xY 1 ^ 2 )" 2 -^)^ (20) 

where Ni are natural numbers. In that case a 'noncommutativity degree' of x % and x J (j > i) will 
be measuring by the following series in Plank constant 



rp3 nt*^ nr-^ /y*3 



+ J>* i<j, :f>:eHM) (21) 



fe=l 



where : f l k J : are sums of ordered monomials (|20|) with arbitrary number of multipliers 

:/?':= E fr Nl - Nn --M Nl ... Nn : (22) 

Afi,...,V n >0;n>0 

and ff?' 1 "' Nn are complex coefficients. 

To introduce dynamics in the algebra of quantum observables JF(M) assume that evolution 
equations for x % have the same form as for the commutative case 

^ = m) (23) 
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but with a 'quantum' derivation X which includes quantum corrections. This derivation is defined 
by its action on the basic elements 



X{&) = ft* :<£:, : g{ :G ^(M). (24) 



k=l 



In the equation (1241) : X* : means the ordered functions X l {x) from the equation (fl6i) where all x 1 
are replaced by x l and : ^ : can be defined as 

:£'•= E 9i Nl - Nn -M Nl ,., Nn :, ^"GC. (25) 

Ni, ... ,7V„>0;n>0 

The commutation relations ( l2Ti) and the quantum derivation ( 1241) include unknown coefficients 
i— n an( j ^i, l— n should be defined. To obtain them we use time independence of the 

relations (|2T|) . In view of the dynamical equation (1231) it means that the following relations hold 

X^'^ - x l xP - E ^ : fk = 0. *<i> i,j = l,...,n. (26) 
fc=i 

Using the properties of a derivation and the formulas (12 ip and (j24j) one can obtain recurrent 
relations connecting unknown coefficients fV> Nl - Nn anc l g^ Nl "' n . The additional conditions for 
receiving these coefficients gives us invariance of the commutation relations (1211) under involution 
and Hermitian property of basics elements inherited from (|15p 

(x*)* = x*. (27) 

Taking into account the explicit expression for quantum derivation (1241) and definition of 
quantum integrals of motion 

^ = Ml) = (28) 
one can obtain these integrals (if they exist) in the form 



■ J k ■- 2^ Nl! ... :Nn > -n>oh ■ IV1 N 1 ...N n ; h G lL - 



(29) 



In the above formula : Iq : is the ordered J where all are replaced by x l . 

Thus, we have described the proposed quantization algorithm which, in principle, gives us 
an algebra of quantum observables with commutation relations ( I2ip among its basis elements, 
a quantum derivation (1241 with related quantum equations of motion (123]) . and a set of quan- 
tum integrals of motion (129]) . Now, let us turn to the examples which illustrate the presented 
quantization method and allow us to clarify some its details. 
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3 Examples 



In this section we apply the proposed method to quantization of some low dimensional classical 
systems. These systems were chosen as examples in spite of the fact that most of them are well 
known. This is because each of the systems has interesting specific features and it is important 
to demonstrate as the method deal with them. In what follows we omit the formal intermediate 
step (transition to commutative algebras) and simple calculations. 



3.1 Charged spinless particle in a magnetic field 

This example demonstrate that one can easily obtain the standard formula by the presented 
algorithm using only a measurable magnetic field without references to a vector - potential (see, 
for instance, [H], [12]). 

A charged spinless particle motion in a uniform magnetic field describes by the Newtonian 
equation 

dv o 

m— = -vxB (30) 
dt c ( J 

where q, m and v are particle's charge, mass and velocity respectively, c is the speed of light, B 

is a magnetic field. For a plane motion (when v is orthogonal to B) the equation ( 1301) gives 



dVx 1 B y ( \ y Q B v i \ fai\ 
-rr = — v y = X m (v x ), — = v x = X m (v y ) (31) 

dt mc dt mc 



where X m is the related vector field. These equations up to notations and interpretation of the 
dynamical variables coincide with the harmonic oscillator equations. So it is naturally expect 
that our quantization procedure will lead to the quantum oscillator algebra. In the case under 
consideration two velocity components v x , v y play the role of the main dynamical variables. Any 
smooth function of the particle velocity v can be considered as an integral of motion. We choose 
here the kinetic energy 

mv 2 r, o 
v 2 = v 2 x + v 2 y . 

Turning to the quantization, let an ordered form of quantum monomials ( l20i) be v^Vy 2 . Then 
using the consistency condition ( |26l) we obtain the following simplest commutation relation 

v y v x = v x v y + ihk (32) 
where k is an arbitrary real number. The quantum derivation 

X m {V x ) = Vy, X m {V y ) = V x 

mc mc 
and the quantum kinetic energy (that is an integral of motion) 



remain without quantum corrections. Let us stress here that the condition (126]) does not give an 
unique solution. In the considered case instead a real number k can be stood an arbitrary function 
of H . To fix this function it is necessary an additional condition. Let the quantum derivation X m 
of the algebra J r (v x ,v y ) be an inner one 

X m (.) = ^[.,H], HeF(v x ,v y ). (33) 
This condition means that the quantum equations of motion have the Heisenberg form 

dv x 1 r , a, dv y 1 ~ 

and uniquely defines the constant k = — 2 -j. Taking into account this value of k we obtain from 
( 1321) the well known quantization rule for the velocity components 

[v x ,Vy] = ih-^-z-. (34) 
m 2 c 

3.2 Euler top 

The example in this section shows that using our procedure one can obtain commutation relations 
of the angular momentum type directly from the classical equations without reference to their 
canonical representation by position and momentum operators. 

The evolution of the angular momentum of a classical anisotropic top with a fixed point at its 
center of inertia is defined by the Euler equations 

= a\L 2 L ?> = X EU (Li) 

= a 2 L x L z = X EU (L 2 ) (35) 
= o, 3 L 1 L 2 = X EU (L 3 ). 

where Li,L 2 , L% are components of the angular momentum vector and real parameters a±, a 2 , 0,3 
depend on principal moments of inertia h, I 2 , 13 by the formulas [10] 



a 1 



h' 1 - h~\ a 2 = /r 1 - v\ % = h' 1 - /r 1 . 



It is convenient regard a 1: a 2 ,a 3 as components of a three-dimensional vector a = {ai,a 2 ,a^). 
Quadratic integrals of motion form a family 

^ = i(6 1 L2 + 6 2 L2 + 6 3 L2) (36) 
where the parameters bi,b 2 , 63 satisfy the relation 

a • b = ai&i + a 2 b 2 + 0363 = (37) 
which means the orthogonality of the vectors a and b = (pi, b 2 , 63). 
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The quantization algorithm gives in this case the following simplest commutation relations 
among the quantum dynamical variables 

L 2 Li = L X L 2 - ihf 3 L 3 , L 3 Li = L1L3 + ihf 2 L 2 , L 3 L 2 = L 2 L 3 - ihfxLx (38) 
and the quantum derivation with the following corrections 

Xeu(Li) = a x L 2 L 3 - 

X EU {L 2 ) = a 2 L l L 3 +'^L 2 (39) 

Xeu{L 3 ) = a 3 LiL 2 — lM ^ A L 3 . 
The real parameters f = f 2 , f 3 ) as well satisfy the orthogonality condition 

a • f = 0. (40) 

Note that the Euler quantum dynamical equations in view of ( 1391) can be rewritten in the sym- 
metrized form 

^Hr = + L 3 L 2 ) 

^ = f(L 1 L 3 + L 3 L 1 ) (41) 

= S 2 (L\L 2 + L 2 Li). 
As in the classical case the family of the quantum observables 

H = i(6 x L? + b 2 L\ + b 3 Ll) (42) 



are the family of integrals of motion if the condition fl37l) holds. If we choose b such that the 
following relation satisfy 

b x f = a (43) 

where x denote a vector product in M 3 , then related integral of motion from the family (|42p 
generates the quantum equations of motion (1411) in the Heisenberg form. We can choose from the 
family (j4"2"j) two independent integrals of motion H, H' which are defined by the vectors b, b' and 
satisfy fl37|) . These integrals are commutative when the following condition is fulfilled 



f ■ (b x b') = 0. (44) 

In this case it is not difficult exercise to obtain Fock space of common states for commuting H and 
H' using, for example, Schwinger oscillator representation of angular momentum type operators. 
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3.3 Fourth order Pais - Uhlenbeck oscillator 

The example in this section shows that our procedure allows to quantize a system in a reasonable 
way whereas its standard quantization leads to the Hilbert space containing negative norm states 
[13], [H]. The Pais - Uhlenbeck oscillator was introduced in the paper [13] as a toy model to study 
field theories with higher derivative terms. The details of the following calculations see in [15j. 
The considered oscillator is defined by the fourth order equation 

^ + H 2 + ul)^ + u 2 x u 2 2 x = 0. (45) 

where Ui,^ uj 2 are positive frequencies. It can be written in the form of a system of first order 
equations 

^ = x 2 = Xp U {x 1 ), ^ = x 3 = X PU (x 2 ) 

(46) 

where x\ = x is the coordinate in the original space, x 2 is the velocity and so on. This system has 
two independent integrals of motion 

Hi = -\J\uj\x\ + \{uj\ + ujI)xI + \x\ + uj\uj\x x x z , 

(47) 

H 2 = \uJ\ijj\x\ + \{ul + u\)x\ - \x\ + x 2 x 4 , 

Our quantization algorithm gives for quantum basis observables the following nonzero commuta- 
tion relations 

[x 1 ,x 2 } = ihf, [x 1 ,x 4 \ = ihg, 

(48) 

[x 2 , x 3 ] = -ihg, [x 3 , x 4 ] = -\hJ{ijjlj - ih(uf + u>l)g 

where /, g are arbitrary parameters. The dynamical equations and integrals of motion have no 
quantum corrections. Put a Hamiltonian in the form 

H = (fujj + g)-\ful + g)-\-fH x + gH 2 ) 

where Hi,H 2 are Hi,H 2 with all Xj are changed by x.i. Then the quantum dynamical equations 
can be written in the Heisenberg form. The appropriate choice of the parameters /, g leads to 
the quantization of the Pais - Uhlenbeck oscillator as the usual anisotropic oscillator with positive 
norm states |15|. 



3.4 Nonlinear oscillator with coordinate - dependent mass 

The example in this section shows that our quantization algorithm allows to quantize nonlinear 
systems. Here we consider an oscillator with coordinate - dependent mass (see [16] and references 
therein). The evolution of this oscillator are defined by the following nonlinear equation 

^l + Xx 2 )^-XxA 2 + u 2 x = 0, A>0 (49) 
at 2 at 
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and it has the non-polynomial integral of motion 



1 y 2 + uj 2 x 2 

2 1 + Ax 

The equation ( 1491) is equivalent to the system 



H = ^-r^r- (50) 



dx dy Xx 2 uj 2 x 

dt =y ' dt = l + Xx 2V ~ 1 + Xx 2 ' ( } 

Using the quantization procedure we obtain the following simplest commutation relation 

[y, x] = -ih(l + Xx 2 ) (52) 
quantum dynamical equations with quantum corrections 

da; A dy Xx ~ oj 2 x 2 2X 2 h 2 x 

Tt =y ' Tt = TTx^ y ~ TTxx-i + lhX{l ~ TTx^ )y + TTx^ (53) 

and corrected quantum integral of motion 

H = l^T-^ 2 + -V) + + h 2 X(-±— - \). (54) 

21 + Xx z 1 + Xx z 1 + Xx z 2 

It is easy to check that H considered as a Hamiltonian allows to rewrite the equations (1531) in the 
Heisenberg form if one take into account the following formula 

1 2ihXx 
^' 1 + Xx 2 ^ = 1 + Xx 2 ' 
This formula can be proof using the integral representation 

i r°° 



1 + Xx 2 

and well known decomposition 

e A Be- A = B + [A,B] + hA, [A, B]] + ... . 

The representation of the obtained algebra of quantum observables with commutation relation 
(1521) see in M- 



We conclude this work by some short remarks. The adduced examples were shown that the 
proposed quantization method is effective and effortlessly gives results in the simple cases. It 
suitable for multi-Hamiltonian systems (the examples 2) and 3)). For Euler top the method gives 
the commutation relations of angular momentum type (|38|) without usual handling of classical 
expressions. In the case of Pais - Uhlenbeck oscillator our algorithm naturally leads to the positive 
norm Fock space of states. The example 4) was shown that the procedure works effectively for 
nonlinear systems as well. On the other hand, it follows from the examples that to obtain the 
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unique result of quantization it is often necessary an additional condition besides the ones were 
declared in Section 2. In the considered cases it was the condition of Heisenberg form of quantum 
dynamical equations. 

The presented examples do not cover most of interesting systems such as dissipative systems, 
system in a curved space-time, etc. We hope that quantization of these systems by the proposed 
method and more precise mathematical formulation of the method itself will be carried out in the 
future. 
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